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Bounds for Hamiltonians with arbitrary kinetic parts 
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A method is presented to compute approximate solutions for eigenequations in quantum mechanics with 
an arbitrary kinetic part. In some cases, the approximate eigenvalues can be analytically determined and 
they can be lower or upper bounds. A semiclassical interpretation of the generic formula obtained for the 
eigenvalues supports a new definition of the effective particle mass used in solid state physics. An analytical 
toy model with a Gaussian dependence in the momentum is studied in order to check the validity of the 
method. 
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I. INTRODUCTION 

The envelope theoryir— is a powerful method to ob- 
tain approximate solutions, eigenvalues and eigenstates, 
of eigenequations in quantum mechanics. It has been 
rediscovered and developed recently under the name of 
the auxiliary field method 7-12 . Both techniques are 
equivalent^, but they were introduced from completely 
different starting points. Let us assume that the Hamil- 
tonian H can be written as (in the following, we will work 
in natural units h = c = 1) 



H=T(p) + V(r), 



(1) 



with p = \p\ and r = \r\. Such a form is relevant for 
one-body and two-body systems. The basic idea is to 
replace this Hamiltonian H by another one H which is 
solvable, the eigenvalues of H being optimized to be as 
close as possible to those of H . Following the structure 
of the Hamiltonian, the approximate eigenvalues i) can 
be upper or lower bounds, or not to have a variational 
character; ii) can be obtained on a closed- form expres- 
sion, or only numerically computed. In the most favor- 
able case, an analytical bound, the dependence of the 
eigenvalues on parameters of the Hamiltonian and the 
quantum numbers can be determined, giving deep in- 
sights about the structure of the solutions and a reliable 
estimation of the spectrum. Even in the less favorable 
situation, a non-variational numerical approximation, it 
is possible to check easily and rapidly more elaborate nu- 
merical computations. 

At the origin, the method has been developed for 
Schrodinger equations^ 7 -, and afterwards it has been ex- 
tended for the semirelativistic kinematics^. The pur- 
pose of this work is to generalize the technique to arbi- 
trary kinetic operators. This is motivated by the exis- 
tence of non-standard kinetic energies in some physical 
problems, for instance in atomic physics (non-parabolic 
dispersion relation)ii or in hadronic physics (particle 



mass depending on the relative momentum) 15 . Another 
motivation is to support the new definition of the effec- 
tive particle mass proposed in Refs i 14 i 16 ~ — . 

The generic method to compute approximate solutions 
(which could be upper or lower bounds) of a Hamilto- 
nian with a non-standard kinetic part is presented in 
Section [TTJ As this work is a generalization of the one 
described in Refj^ 2 -, the main equations are worked out 
without too many details. In Section IIII1 a semiclassi- 
cal interpretation of the generic equations is given. In 
order to check the validity of the method, an analytical 
toy model with a Gaussian form for the kinetic part is 
solved and the formula obtained is compared with nu- 
merical solutions in Section llVl Concluding remarks are 
given in Section IVl 



II. GENERAL EQUATION 

The envelope theory is generalized here to treat on the 
same footing the potential and kinetic parts. The idea is 
to replace the Hamiltonian (TTJ) by the following one 



H = T + V 1 



with 



and 



T = 



El 

2v 



■T{J{v))- 



JJyf 
2v 



V = pP(r)+V(I(p))-pP(I(p)). 



(2) 



(3) 



(4) 



v and p are two real parameters, and we assume that the 
following functions are well-defined 



I{x)=K~ x {x\ K{x) 



P'{x) ' 



J(x)=L- 1 (x), L(x) = - 



T'(x) 



(5) 



a ' E-mail: claude.semay@umons.ac.be 



The kinematics of H is nonrelativistic, P(x) is an aux- 
iliary potential, and a prime denotes the derivative. An 



2 



eigenvalue of this Hamiltonian is given by 



E{v,p)=T(J(u)y 



Jiy) 

2v 



■ + V{I{p))-pP{I{p)) + e{v,p), 

(6) 



where e(v, p) is an eigenvalue of the nonrelativistic Hamil- 
tonian 



H NR = ^+pP(r). 



(7) 



Kinetic and potential parts are treated in a similar way, 
but with x 2 which is the counterpart of P(x) and 1/(2 v) 
which is the counterpart of p. The approximation for 
an eigenvalue of the genuine Hamiltonian is given by 
E(vo,Pa) for which 



dE(u,p) 



dE{v,p) 



dp 



= 0. 



(8) 



Within these conditions, T is tangent to T and V is tan- 
gent to V. The comparison theore m 19 i 20 implies that, if 
T > T and V > V for all values of the arguments, then 
the eigenvalues of H are upper bounds of the correspond- 
ing eigenvalues of H . Reciprocally, if T < T and V < V 
for all values of the arguments, then the eigenvalues of 
H are lower bounds. In other cases, no guarantee exists 
about the variational character of the approximations. A 
simple criteria to determine if a bound exists is given in 
Refsi£~— for the potential, but it can be also used for the 
kinetic part. Let us define two functions h and g such 
that 



T(x) = h(x 2 ) and V{x) = g(P(x)). 



(9) 



If h"(x) and g"{x) are both concave (convex) functions, 
E(v ,p ) is an upper (lower) bound of the genuine eigen- 
value. If T(p) oc p 2 (V(r) oc P(r)), the variational char- 
acter is solely ruled by the convexity of g(x) (h(x)). 

Interesting results can be obtained if the auxiliary po- 
tential is a power law, 

P(r) = sgn(A) r A with ^ A > -2. (10) 

Within this condition, p and v are always positive quan- 
tities, and e(y, p) can be written under the formal 7 - 



eO, P) 



2X A/(A+2) 



(11) 



where Q is a global quantum number. The method is 
particularly interesting if Q is exactly known. This is the 
case for the Coulomb interaction (A = — 1, Q = n + l + l) 
and the harmonic potential (A = 2, Q = 2n + I + 3/2). 
If A = 1, Q is known only for I = states and is equal 
to 2(— a n /3) 2 / 3 , where a n is the (n + l)th zero of the 
Airy function Ai. For arbitrary values of A, simple and 
good analytical approximations can be found in Ref. 7 . 
But, if Q is not computed with a sufficient accuracy, the 
variational character of a bound cannot be guaranteed. 



After some algebra, constraints ((SJ, with P(r) given 
by PHD ; reduce to 



(|A|Q A )^2(p ^o) TT5 =Po, 



\M Q 



: {po vo) 



|A|r, 



(12) 
(13) 



where r = I(po) an d Po = J(vo)- From (|T2|) and (fT3|l . 
we can deduce that ropo — Q- Taking into account this 
result, plus the relations p — K(r ) = V'(ro)/(|A| r^ 1 ) 
and v = L(po) = p /T'(p ), (fT2|) and (fl~3l) can be writ- 
ten in the more compact form (|16p . After some algebra, 
the approximate eigenvalue E(v ,p ) given by ^ with 
the parameterization (jlip can be greatly simplified into 
(jT4j) . Finally, the approximate solution is given by the 
following set of equations 



E = T(po) + V(r ), 
Q 

ro ' 

PoT'(po) =r V'(r ). 



Po 



(14) 
(15) 
(16) 



The parameter ro can then be interpreted as a mean dis- 
tance between the particles and po as a mean momentum 
per particle. Both parameters depend on the quantum 
numbers via Q. The value E can be a (upper or lower) 
bound following the convexity of functions h(x) and g{x), 
as explained above. Let us note that the only trace of 
the auxiliary potential is contained in the value of Q, and 
that (1161) is the translation into the variables ro and po 
of the generalized virial theorem^!. The system (|14|) -(fl6 |) 
is similar to the systems (3.2)-(3.4) in Ref.— and (15)- 
(17) in Refi^, but here the form of the kinetic part is 
arbitrary. 

Let us note | f , p ) an eigenstate of the nonrelativistic 
Hamiltonian i?NR given by ([7]) with v = vq = pa/T'{po) 
and p = pa = V^'(r )/(|A| r^ 1 ). Such a state is an ap- 
proximation of the corresponding eigenstate of H 9 . Using 
the Hellmann-Feynman theorem 2 ^ as in ReL 7 -, it can be 
shown that 



(v ,p \p 2 \u ,po) =pl, 
(vo,po\r X \is ,Pa) = Tq. 



(17) 
(18) 



This confirms the interpretation of parameters ro and pq 
as mean values. 



III. SEMICLASSICAL INTERPRETATION 

In Refs . 14 ' 16 " — , an effective particle mass m eff is de- 
fined by the relation p = m cS v, where v is the group 
velocity of the associated wave packet. It follows then 
that 



off (dT 



(19) 



Using this definition, a semiclassical interpretation of the 
system (fT4]) - P^|) is possible. If Ti(pi) is the kinetic energy 
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m| ff 



Figure 1. Classical circular motion of two particles in their 
center of mass frame. 



for the ith particle, the kinetic energy for two particles 
can be written Ti(po) + ?2(po) = T(po), where po is the 
module of the common momentum in the center of mass 
frame. With the effective mass (1191). we have 



off 



Pi) 



(20) 



and the speed of the ith particle is given by Uj — T/(p )- 
As in Ref. 12 , let us assume a classical circular motion for 
the two particles (see Fig. Q} . The centripetal force F{ 
acting on the zth particle is given by 



F, 



-cff ' 



— = Pa- 
ri 



TRpo) 



(21) 



If ro is the distance between the two particles, the rigid 
rotation constraints, ro = r\ + r 2 and v\/ri — vi/ri, 
imply that 



ro 



T'(Po) 



(22) 



If the force acting on the zth particle comes from the 
potential V(r) generated by j, then F\ = F 2 = V'(r ). 
(|2Tj) and can be recast into the form (flU)) , and it 
is obvious than (TT4)) gives the mass of the system. A 
semiclassical quantification of the total orbital angular 
momentum gives ro po = 1 + 1/2, and we obtain a system 
very similar to (fT^ - pT)]) . This supports P^)l as a good 
definition for the effective mass. 



IV. A TOY MODEL 

In this section, we solve a simple toy model in order 
to check the relevance of the method. Let us consider 
the following Hamiltonian with a Gaussian form for the 
kinetic part 



where the parameter m plays the role of a mass. Indeed, 
for high values of m (3> a 1 / 3 ), H reduces to a harmonic 
oscillator Hamiltonian 



H 



2 777 



a r 



0(4 



(24) 



The parameter a = 1 or 2 is the number of particles 
(arbitrary positive value of a can also be considered 
to study duality relations between different many-body 
systems^). Such a Hamiltonian is not realistic but it 
has been chosen because i) it admits an analytical lower 
bound (see below); ii) it reduces to the well known case 
(|24p in a well-defined limit; Hi) accurate numerical solu- 
tions are not easy to obtain (see below). It is more conve- 
nient to work with the conjugate dimensionless variables 
x = yja/(a m ) r and q = -J a m/ap for the dimension- 
less Hamiltonian = H/ (a m) given by 



H A = cxp (k q 2 ) + x 2 with k = 



2 a 77i 3 



(25) 



The corresponding eigenvalues are noted e = E/(am). 

e(k) 
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Figure 2. Eigenvalues e of (|25[l as a function of k. By in- 
creasing energy, the eigenvalues corresponds to (n, I) = (0, 0), 
(0,1) and (1,0). Dots: accurate numerical solutions; Solid 
black lines: lower bound (|26|) with Q = 2n + I + 3/2; 
Dashed grey lines: harmonic oscillator approximation (|27[) 
with Q = 2ra + Z + 3/2. 



Using the set of equations (|I^ -(|T6 j) . the following ap- 
proximate solution can be found 



exp ( 2 W [ Vk 



Q 



1 + 2 Wo { Vk 



(26) 

where W'o(z) is the main branch of the Lambert function 
(also called Omega function or product logarithm)^. If 
k <C 1, we can write e app = £ho + O (k) with 



e HO = l + 2v / fcQ. 



(27) 



H = a 777 exp I ^ } + ar 2 

2 77i/ 



(23) 



This corresponds to the harmonic oscillator approxima- 
tion. A natural choice is to take P(x) = V(x) = x 2 . 
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Then, Q — 2n + 1 + 3/2 and the function h(x) defined by 
© is convex. So, (j2"6]l is a lower bound, whose quality 
is shown on Fig. [2] The numerical solutions of (|25|) have 
been computed with the the three-dimensional Fourier 
grid Hamiltonian method 2 ^— which is particularly well 
suited for this type of Hamiltonian. This numerical pro- 
cedure is equivalent to an expansion in a special basis 27 
and implies the computation of Hamiltonian matrix ele- 
ments. Because of the exponential function in (|25l) . these 
matrix elements can be huge numbers and they must be 
computed with a very high accuracy to obtain stable and 
accurate eigenvalues. One can see that the lower bound 
is quite good and more accurate that the harmonic oscil- 
lator approximation when k increases. 



V. CONCLUDING REMARKS 

The envelope theory (or equivalently the auxiliary field 
method) is a method to compute approximate solutions 
(generally upper or lower bounds) of Hamiltonians of the 
form H = T(p) + V(r), where the kinetic part T is a 
nonrelativistioi»Z or a semirelativistic one&I^. In this pa- 
per, it is shown that the method can be used for arbitrary 
forms of T. The idea is to replace the Hamiltonian H by 
another one H = T + V which is solvable, and with T 
and V respectively tangent to T and V. Provided T 
is a nonrelativistic kinetic operator and V a power-law 
potential, the approximate eigensolutions can be easily 
computed by solving a set of three equations which have 
a natural semiclassical interpretation. Nevertheless, the 
computation is a full quantum one since eigenvalues and 
eigenstates are obtained for a well-defined global quan- 
tum number Q. 

A priori, This method can be applied to a wide vari- 
ety of Hamiltonians relevant for various domains, from 
atomic to hadronic physics. With a good choice of the 
power law, the value of Q is analytically known and the 
eigenvalues obtained can be upper or lower bounds of 
the genuine energies. If numerical approximations can 
be easily computed, closed-form formulae can even be ob- 
tained for some particular Hamiltonians. The toy model 
studied here with a harmonic potential and a Gaussian 
kinetic operator is in the most favorable situation: A 



quite accurate analytical lower bound is obtained giving 
deep insights about the structure of the solutions. 
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